Dynamic Ordering and Transverse Depinning of a Driven Elastic String in a 

Disordered Media 



O 

o 

(N 



in 



> 
in 
o 

O 



o 
o 



X 



C. ReichhardtW and C. J. OlscV 2 ) 
CNLS and Applied Theoretical and Computational Physics Division, 
' 2 - ) Theoretical and Applied Physics Divisions, 
Los Alamos National Laboratory, Los Alamos, NM 87545 
(February 1, 2008) 

We examine the dynamics of an elastic string interacting with quenched disorder driven perpen- 
dicular and parallel to the string. We show that the string is the most disordered at the depinning 
transition but with increasing drive partial ordering is regained. For low drives the noise power is 
high and we observe a l// 2 noise signature crossing over to a white noise character with low power 
at higher drives. For the parallel driven moving string there is a finite transverse critical depinning 
force with the depinning transition occuring by the formation of running kinks. 
PACS: 64.60.Ak, 71.45.Lr, 74.60.Ge 



The dynamics of driven elastic media interacting with 
quenched disorder occurs in a wide variety of physical 
systems which include magnetic domain wall motion jlj , 
nonequilibrium growth models of friction |||, vortex 
lattice motion in superconductors and charge den- 
sity waves Recently an intense interest in driven 
elastic media with quenched disorder has been directed 
at dynamical ordering where at low applied drives the 
system is in a highly disordered pinned state and at a 
critical drive a depining into a highly disordered moving 
state occurs, while at high drives the effect of the disor- 
der is reduced and the system can regain a considerable 
amount of order. Particular systems where this reorder- 
ing behavior has been studied extensively experimentally 
§-£|, theoretically fl^lf], and numerically @|§Q 
include vortices in superconductors as well as charge- 
density wave (CDW) systems. 

In the vortex system if there is no quenched disorder 
the vortices form an ordered crystalline state; however, 
in most samples there are defects which attract and pin 
vortices, disordering the vortex lattice. When a driving 
force is applied, the initial depinning transition can be 
plastic, with certain regions of mobile vortices tearing 
past pinned regions and a large number of defects gener- 
ated so that the system has only a liquid like structure. 
For higher drives the vortices regain order; however, the 
pinning can still affect the structure so the vortices can 
have a partial anisotropic ordering or smectic structure 
|jl2f . An open question is whether the reordering transi- 
tion is a true transition or a crossover. 

An interesting feature of the reordered moving state is 
that a transverse depinning threshold can exist as orig- 
inally proposed in Ref. pi] . Here, although the lattice 
is moving in the longitudinal direction, the effect of the 
pinning is still present in the transverse direction so that 
a finite transverse drive must be applied before the lattice 
will move in the transverse direction. Simulations [ p~4|j3~T| 
and experiments Q on vortex lattices as well as simu- 
lations for friction models || have found evidence for a 
finite transverse depinning threshold for the longitudinal 



moving systems. Other theoretical work contends that 
at finite T no true critical transverse depinning thresh- 
old exists; instead, due to the high pinning barriers, a 
pronounced crossover in the transverse IV curves may 
appear [|l2| . In the CDW context the effect of the quasi- 
particle current flowing perpendicular to the CDW 2k p 
ordering wavevector was also studied theoretically |i"§| l 
and in experiment [fl9f . An intriguing question is whether 
the reordering transitions and the transverse barrier seen 
in the vortex system occur in other types of driven elastic 
media, such as driven interfaces and polymers interact- 
ing with quenched disorder. An ideal way to model these 
systems in 2D is with a driven elastic string interacting 
with quenched disorder. 

Previous numerical studies of driven elastic strings 
with disorder [p0|-p3| employed drives perpendicular to 
the string and focused on the dynamics near the depin- 
ning threshold. In this case critical behavior is expected, 
and was observed in the form of broad distributions of 
avalanche sizes as well as scaling of the velocity near de- 
pinning, v = (/ — f c ) a , as proposed by Fisher [ pi[ . The 
roughness of the string as measured by the roughening 
coefficient exponent was 0.9 to 1.1 at depinning [^0 22 
and 0.5 at high drives fl20|| , indicating that some ordering 
of the string was occurring as the drive increased. 

In this work we investigate the continual evolution of 
the order and noise signatures in driven strings for both 
parallel and perpendicular drives. We also investigate the 
transverse depinning for the moving string. In the case 
of the perpendicular moving string we do not observe a 
transverse barrier; however, for the parallel driven string 
a transverse depinning threshold is observed which oc- 
curs through the formation of running kinks. For in- 
creasing longitudinal drives the transverse barrier de- 
creases. Physical systems relevant to the perpendicular 
driven string include domain walls and moving interfaces, 
while systems similar to the parallel driven string include 
boundaries between sliding surfaces, polymers that are 
aligned with an applied drive in a random media, or a 
single stream of fluid flowing down a rough surface. 
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FIG. 1. The string configuration at different drives for a 
system with driving in the x-direction perpendicular to the 
string with a depinning threshold of F dp = 0.09. (a) Right 
above depinning at F d /F dp = 1.1; (b) F d /F dp = 1.38; (c) 
F d /F dp = 1.94; (d) F d /F dp = 4.4. 

We model an overdamped elastic string in (1 + 1) di- 
mensions. The string is composed of discrete particles 
connected by springs. The system has periodic bound- 
ary conditions in the x and y-directions and the string is 
connected periodically in the y-direction. The equation 
of motion for each particle on the string is: 



= Kfs +t P + Fd 



(1) 



Here r is the particle location, 7 = 1 is the damping term, 
k = 5 is the string elastic constant, f s = r»_i — 2rj + rj + i. 
is the spring force from the two neighboring particles, f p 
is the pinning force and is the uniform applied driv- 
ing force. The pinning is modeled as attractive parabolic 
traps scattered randomly through the sample. The pin- 
ning force is f p = - Ef P (/p/ r p)( r * - r i P) ) e ( r p ~ \ Y * ~ 
rj^l), where O is the Heaviside step function, rjjf-' is the 
location of pinning site k, f p = 0.3 is the maximum pin- 
ning force, and r p = 0.25 is the pinning radius. The 
results shown here are for strings containing A" = 500 
to 2000 particles, interacting with N p — 4.6 x 10 5 pin- 
ning sites in samples of size 160 x 160. We consider two 
types of initial conditions: one where the string is put 
down in its unstretched equilibrium position, and a sec- 
ond where the string is annealed at a finite temperature 
by the addition of a thermal kick. We find that when the 
applied drive is increased slowly enough, both methods 
give similar results. We apply a uniform force Fd in the 
x direction for perpendicular driving or the y-direction 
for parallel driving. We increase the driving force in 
increments of 0.001 and spend 30000 simulation steps 
at each increment, measuring the average string velocity 
V x = x • Vi, and V y = ^f* y • Vj. To quantify the 
order in the string we measure the difference AL in the 
length of the string compared to the equilibrium length, 
AL — X^G 1 "* — r i+i| — r 0> where r; is the equilibrium 
length of each string segment. We also measure the noise 
in the string velocity: S(w) = \ J V x (t) e - l27Tult dt\ 2 , and 
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FIG. 2. (a) The string velocity V x versus F d for the system 
in Fig. 1. The string is pinned for F d < 0.09. (b) AL ver- 
sus F d showing that the string becomes disordered as F dp is 
approached from below, and gradually reorders for increasing 
F d > F dp . (c) Noise power So versus F d . 

the noise power in one frequency octave So = S(ui). 

In Fig. l(a-d) we show a series of snapshots of a per- 
pendicular driven string. In Fig. 1(a) at a driving force 
Fd/Fdp = 1.1 very close to the depinning transition 
Fdp = 0.09, the string is highly disordered. Above de- 
pinning, the motion of the string is jerky, occurring in 
pulses or avalanches, with portions of the string immobile 
while other portions are moving as can be seen from the 
instantaneous velocity of individual beads. In Fig. 1(b) 
and (c) for Fd/Fd p — 1.25 and 1.75, the string becomes 
less winding. Finally in Fig. 1(d), for F d /F dp = 4.0 the 
string is almost completely ordered again. 

To quantify the evolution of the driven string from 
Fig. 1, we plot in Fig. 2(a) Fd versus < V x > show- 
ing that the string is pinned for Fd < 0.09. In this case 
the string was initially placed in its equilibrium position. 
In Fig. 2(b) we plot AL which shows that the string is 
aligned for Fd — but gradually disorders as F^ p is ap- 
proached from below. AL reaches a maximum at depin- 
ning. For increasing F d > F dp , AL gradually decreases as 
the string orders. In Fig. 2(c) the reordering can also be 
seen in the noise power So which is maximum at depin- 
ning and again decreases for increasing Fd- The decay of 
5*0 fits well to an exponential form where for Fd > Fd p , 
So(Fd) oc exp(— (Fd — Fdp))- We have also conducted 
simulations where we decrease Fd after a ramp up, and 
we do not observe any hysteresis in the ordering AL or 
the noise power So. 

In Fig. 3 we show noise spectra for the same system. 
Right above depinning, at Fd/Fd p = 1.1, S(u>) shows a 
l// 2 characteristic, as seen in Fig. 3(a). For increasing 
drives, the string begins to move continuously and re- 
order. As a result, a characteristic time scale forms, the 
l// 2 shape is lost and the higher frequencies are reduced 
in power, as shown in Fig. 3(b) for Fd/Fd p = 4.4. Noise 
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FIG. 3. (a) Noise spectra for the system in Fig. 1 

for Fd/Fdp — 1.1 showing a l// 2 signal. (b) S(u>) for 
Fd/Fdp = 4.4 showing a characteristic time scale. 

spectra for vortices have shown similar features in both 
experiments p5fl and simulations [ fL6| . In the vortex sim- 
ulations the l// 2 high noise region corresponds to the 
disordered plastic flow state and the low noise power re- 
gion corresponds to the reordered elastic flow state. 

We have repeated the same type of simulations shown 
in Figs. 1-3 for the parallel driven string and observe 
the same general features. For the parallel driven string, 
the depinning threshold is lower and we observe a clear 
washboard signal for the high drive case. We have carried 
out a series of simulations for varied f p , k, and N p . In 
each case, the depinning threshold changes; however, the 
same general features of the reordering are observed. 

To test the predictions in Ref. |ll|] of a critical trans- 
verse threshold for the driven string, we have conducted 
a series of simulations where we fix the driving force at 
a constant value after the string is in motion. We then 
slowly apply an additional driving force in the direction 
transverse to the initial longitudinal motion. For the per- 
pendicular driven string the applied transverse drive is 
along the string and for the parallel driven string the 
applied transverse drive is perpendicular to the string. 
For the perpendicular driven string we do not observe a 
transverse barrier for the moving string. For a parallel 
driven string, in Fig. 4(a) we show the transverse driving 
force versus the transverse string velocity indicating that 
for the case of the parallel moving string there is a finite 
transverse pinning threshold. We note that for the mov- 
ing vortex systems a finite transverse barrier occurs only 
when the vortices are moving in well defined channels 
and the transverse drive is perpendicular to these 
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FIG. 4. (a) The transverse string velocity V versus the 
transverse drive Fj for a parallel driven string (lower curve). 
The parallel string drive is Fd/Fd p = 1.22. A finite transverse 
depinning threshold is visible. A plot of V for a system with 
no quenched disorder is also shown (dashed curve), (b) The 
transverse depinning threshold Fj p versus longitudinal drive 
F d . The solid line indicates a fit to Fj p oc (F d — F dp )~ 2/:i 

channels. The parallel moving string may be viewed as 
a single moving channel. In Fig. 4(b) we show that 
the transverse barrier decreases with increasing longi- 
tudinal velocity. In the vortex system as the longitu- 
dinal drive increases the channels straighten out and 
the transverse barrier decreases. In the parallel moving 
string the decrease in the transverse pinning threshold 
also corresponds to the string reordering for the higher 
longitudinal drives. Above the parallel depinning tran- 
sition, Fd > Fd p , the data follows a power law form, 
Fj p oc (F d - F dp )~ 2 / 3 . In a CDW system an exponential 
decay would be expected but for the driven string 
system a power law decay is predicted instead Q. 

We have also investigated the dynamics of the trans- 
verse depinning for the parallel driven string, as indi- 
cated in Fig. 5, which shows a series of snapshots of the 
string just above the transverse depinning threshold. The 
transverse depinning transition is very distinct from the 
longitudinal depinning of the perpendicular driven string 
shown in Fig. 1. The transversely driven parallel mov- 
ing string does not become more disordered as observed 
at the depinning transition of the strictly perpendicular 
driven string. Instead, the transverse depinning occurs 
through the formation of a running kink that moves in 
the direction of the longitudinal drive as seen in Fig. 5. 
Although one might expect a kink and an anti-kink to 
form, we observed that below the kink the string is po- 
sitioned at a slight angle. The kink moves along the 
string through the periodic boundary conditions. As the 
transverse drive is increased more kinks appear and the 
kinks begin to move more rapidly. Similar kink motion 
has been observed in interfaces in sliding friction systems 
. The moving kinks themselves can show intricate dy- 
namics. We have observed that these kinks can become 
pinned, and that another kink can form elsewhere and 
later collide with the pinned kink to form a larger kink. 

In conclusion we have numerically investigated the dy- 
namics of perpendicular and parallel driven strings in 
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FIG. 5. (a) A series of snapshots of a parallel driven 
string as a function of time with the longitudinal drive in the 
{/-direction and the transverse drive in the ^-direction, for a 
system with Fd — 0.12, just above the transverse depinning 
threshold. The transverse depinning occurs via formation of 
a running kink that moves in the direction of the longitudinal 
drive, (b) The same snapshots offset in the x direction. 

quenched disorder. We find that the string is most disor- 
dered at the depinning transition but recovers order for 
increasing drive as measured by the difference between 
the length of the string and the equilibrium length. This 
reordering is observable in the noise characteristic. Broad 
l// 2 noise is observed at depinning, with a gradual reduc- 
tion in the noise power for increasing drive and a shift to 
a whiter noise characteristic. For the parallel driven mov- 
ing string we find that there is a finite transverse critical 
depinning threshold Fj p , indicating that the theory of 
Rcf. pd| | regarding the transverse barrier in a very differ- 
ent system still appears to apply to this system. We see a 
decrease of Fj p with increasing longitudinal drive, as pre- 
dicted in Ref. Q for CDWs, and find that the decrease 
follows a power- law form as predicted |^6| . The trans- 
verse depinning occurs by the formation of longitudinal 
kinks. We do not observe a transverse depinning thresh- 
old for the perpendicular driven string. Possible experi- 
mental systems in which to investigate the predications 
in this work include noise measurements and domain wall 
imaging for magnetic domains driven at increasing mag- 
netic field sweep rates, as well as varied moving contact 
line speeds over a rough surface for increased drive. Our 
prediction for the reordering and the transverse barrier 
for the parallel driven string may be observable in poly- 
mers aligned with the drive in disordered media or a sin- 
gle stream of liquid flowing over a rough surface at various 
longitudinal and transverse tilt angles. 
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